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Personnel 

Roger Chevalier and Craig Sarazin (both Professors of Astronomy at the University of 
Virginia) were the Principal Investigators. Co-Investigators were Steven Balbus (Associate 
Professor), Robert O’Connell (Professor), and Joel Dregman (Scientist at NRAO in Char- 
lottesville). At a part of the development of the Virginia Institute for Theoretical Astron- 
omy (VITA), John Ilawley was hired as an Assistant Professor. Post-doctoral fellows who 
were fully supported by the NASA Theory Program Grant included Peter Decker (I h.D. 
University of Colorado), Kamizierz J. Borkowski (Ph.D. University of Colorado), Robert 
Emmering (Ph.D. University of Virginia), Alexander Kashlinsky (Ph.D. Cambridge Univer- 
sity), Stephen Murray (Ph.D. University of California, Santa Cruz), and Noam Soker (Ph.D. 
Technion, Israel). In addition, Eric Myra (Ph.D. University of Pennsylvania, post-doc 'at 
Stony Brook) was hired jointly by VITA and the Commonwealth Center for Nuclear and 
Particle Physics at the University of Virginia. 


Graduate Students 

Graduate students who were supported in part by the grant included Geary Albright, 
Bill Dalton, Robert Emmering, Prudence Foster, Chris Graney, John Houck, Eric Lufkin, 
John Miller, Phil Plait, and Michael Wise. Foster, Houck, and Wise completed their Masters 
Degree last year, and Emmering completed his Ph.D. 


Undergraduate Students 

Greg Ashe and Jay Doisscau, undergraduate astrophysics majors, did his senior theses 
on research supported by this grant. In both cases, the research lead to published papers. 
Ashe was supported during part of the summer following his senior year; during this time, he 
wrote up the work for publication. Both were also supported by the NRAO REU program. 
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Research Results 


Supernovae and Supernova Remnants 

K. Borkowski (post doc on the grant), J. Blondin (post doc), and Sarazin developed 
a detailed model for Kepler’s supernova remnant (SNR). Observations of this SNR have 
revealed a strong interaction with the surrounding circumstellar medium, winch was studied 
through both analytical and numerical (with a 2-D hydro code) calculations. 

Borkowski and M. Shull (Univ. of Colorado) continue to investigate effects of electron 
thermal conduction on the structure of radiative interstellar shock waves. Their primary 
goal is to explain the observed line emission from metal-rich ejecta in supernova remnants, 
incorporating up-to-date atomic data. 

Chevalier and R. Emmering (a post-doc on the NASA grant) continued their investigation 
of light echoes around SN1987A. Analyzing infrared and scattered circumstellar fight echoes, 
they used early observations to set limits on the mass of circumstellar dust. Early in 198. , a 
light echo was observed that was sufficiently close to the supernova to suggest a circumstellar 
(i.e. presupernova mass loss) origin. Chevalier and Emmering found that the echo properties 
were consistent with the shell created by the deceleration of the red supergiant wind by 
the surrounding pressure if that pressure is low. More recent observations suggest a more 
complicated morphology and it is plausible that winds from nearby early type stats play 
a role. In particular, the wind from star 2 may have swept out the slow wind from the 
supernova progenitor in our direction. 

Chevalier and C. Fransson (Stockholm Observatory) completed work on the emission 
from heavy element gas ejected in the supernova explosion of massive stars. 1 bey assumed 
that a radioactive energy source is present and calculated the detailed heating and ionization 
of the gas. At times later than 200 days, the emission is dominated by neutral and singly 
ionized lines, the strongest due to [01], [Call], Mgl], [CIl, and [Sil]. 1 he relative and absolute 
line strengths depend on both the column densities of the different abundance zones and the 
relative abundances within each zone. Reliable modeling requires a knowledge of the density 
structure and the degree of element mixing; emission line profiles are useful for determining 
these. Data on the late emission from Type lb supernovae suggests that these are the 
explosions of massive star cores and are not exploding white dwarfs. 

Chevalier investigated the evolution of supernova remnants in the very high pressure 
environment of a starburst galaxy. 

Chevalier modeled the recently observed increase in radio flux from SN 19S7A as the 
result of the interaction of the supernova shock wave with the termination shock of the wind 
from the blue supergiant progenitor star. 


Chevalier, Blondin, and Emmering completed the first part of a study of hydrodynamic 
instabilities in supernova remnants. 

Chevalier and E. Liang (Lawrence Livermore Natl. Lab.) examined the interaction of 
a massive star supernova with a circumstellar shell resulting from previous mass loss. Self- 
similar solutions can describe several aspects of the flow, 'Die supernova remnant Cas A is 
probably a good example of compact shell interactions. 'Die model predicts that the bright 
X-ray and radio ring is shocked dense wind gas, and hot., shocked supernova gas is inside of 
the ring. Compact radio sources observed in NGC 1119 and in the starburst nuclei of M82 
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and NGC 253 may also be compact shell interaction. The model predicts slow expansion 
of the radio sources and the appearance of new sources with diameters comparable to those 
of the present sources. This interpretation requires that the starburst nuclei have a large 
population of massive star supernova progenitors in a blue evolutionary phase, lhe model 
can also be applied to the interaction of SN1987A with its circuinstellar shell, expected to 
occur within decades. 

Chevalier and N. Soker (a post-doc on the grant) calculated models for the expansion of 
a supernova envelope like that thought to be present in SN1987A. They showed analytically 
that the expanding envelope can be divided into an inner region with a flat density profile and 
an outer sleep power law density profile. For SN1987A, the bend in the profile occurs at a 
velocity of about 4000 km s -1 . Polarization and imaging observations show that the envelope 
expansion of SN1987A is asymmetric, so Chevalier and Soker modeled on asymmetric initial 
density distribution or energy deposition. The analytic and numerical calculations show that 
the final degree of asymmetry is comparable in magnitude to the initial asymmetry. 1 he 
models also show that the sense of density asymmetry changes across the bend m the density 
profile. This accounts for the different polarization angles observed in the Ho and Ca 11 lines 
by Cropper et al. in June 1987. The cause for the asymmetry is uncertain. If the progenitor 
envelope was rotationally flattened, angular momentum transfer from a binary companion 
during the stellar evolution is probably necessary. An asymmetric initial explosion is another 
possibility. 

Emmering and Chevalier used the observed light curve and temperature evolution, as 
well as the theoretically predicted initial ultraviolet burst of SN 1987A to construct models 
for the infrared (1R) echo and scattered light echo arising from circuinstellar dust formed 
in the wind from the supernova progenitor’s red supergiant phase. The echo models were 
constructed by including the detailed emission and scattering of dust grains, and a spectrum 
of grain sizes were considered. The comparison of the predictions of the lit echo models 
with the available 111 observations led to the conclusion that an III echo has probably not 
been observed as of late 19S7 and it was possible to set limits on the amount of circuinstellar 
dust. The methods which were developed to study the effects of circuinstellar dust around 
SN 1987A should applicable to observations of other supernovae. 

Emmering and Chevalier also investigated the intrinsic luminosity and initial period of 
pulsars by using Monte Carlo techniques to generate a Galactic sample of pulsars having par- 
ticular distributions of initial periods, magnetic dipole moments, and positions. Hy choosing 
a relationship between the radio luminosity, the period, and the period derivative, it was 
determined which pulsars of the generated sample could be detected in one of the major pul- 
sar surveys. The properties of the detectable pulsars were compared with the properties ol 
the observed pulsars, and the initial distributions and luminosity relation of the model were 
adjusted until a satisfactory agreement was obtained. A crucial point is that the observed 
relation between luminosity and period is affected by selection affects and does not apply to 
the actual pulsar population. It was found that the best fils were obtained with models in 
which the pulsars are injected with relatively long periods, about O.o sec. 

Observations of the supernova 1987A show the presence of asymmetry in the expanding 
envelope. Hawley and G. Lindahl (graduate student) investigated supernova models that 
have asymmetric initial density distributions or asymmetric explosions through simulations 
of nonspherical explosions and point explosions in rotationally flattened envelopes. 

J. Houck (a graduate student on the grant) and Chevalier continued work on the hydrody- 
namics of optically thick, spherically symmetric neutron star accretion in supernovae. 1 hey 
obtained detailed solutions for the steady state accretion flow including Newtonian gravity, 
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a detailed neutrino cooling function and pressure contributions from radiation, electron- 
positron pairs and baryons. Radiation pressure dominates in most of the post shock icgion, 
although pairs become important once the temperature rises above T » 5 x 10 9 K. The post 
shock flow is nearly adiabatic except for a thin zone near the neutron star surface where 
neutrino cooling is important. They also carried out a linear stability analysis of the ac- 
cretion flow. Unfortunately, difficulties in solving the resulting system of still differential 
equations severely limited the range of accretion rates for which solutions could be obtained. 
The solutions that were obtained show that the flow is stable for accretion rates typical of 
very early times (M ~ 10 6 Af©/yr) and that the flow becomes less stable as the accretion 
rate decreases. Hydrodynamic modeling was inconclusive due to numerical difficulties. 

M. Wise (a graduate student on the grant) and Sarazin calculated the effect of charge 
transfer reactions on the X-ray line emission of supernova remnants. Recent laboratory 
plasma physics experiments have shown that charge transfer rales between neutral hydrogen 
and very highly charged ions are extremely rapid. 1 he charge transfer reactions II + X » 
H+ occur primarily to levels with n « Z/ 2, for which the reactions are nearly 

resonant. These reactions are important in fusion devices, because the gas may be cool 
enough near the walls to allow a significant amount of neutral hydrogen. In addition, many 
fusion devices utilize heating by neutral beams, which results in neutrals being injected 
directly into the hot plasma. In calculating the X-ray emission from supernova remnants, it 
is conventional to include only ionization, recombination, and excitation processes associated 
with collisions between ions and electrons. If a small proportion of the hydrogen behind the 
shock remains neutral, then charge transfer reactions might afTect the ionization state and 
X-ray line emission of the gas. The presence of neutral hydrogen in the shock requires that 
the hydrogen entering the shock not be fully ionized. For shocks associated with ljpe 1 
supernovae or O-star stellar winds, it is likely that the surrounding gas is photo-ionized 
by the present or pre-supernova early-type stars. However, for Type I supernova or I ype 
II supernova in dense regions, the surrounding gas may well be neutral. Lvidcnce that 
this is indeed the case comes directly from the observation of broad and narrow Rainier 
line filaments in a number of supernova remnants, which are produced by charge transfer 
between rapidly moving protons behind the shock and neutral hydrogen atoms entering the 
shock. Wise and Sarazin find that charge exchange can affect the X-ray spectra of supernova 
remnants under some circumstances, but that the effects are not generally very significant. 


Interstellar Medium 

Balbus and Borkowski concluded a study of the time-dependent non-equilibrium ion- 
ization structure and evolution of radiative magnetized plane-parallel thermal conduction 
fronts in the interstellar medium. I he planar cloud geometry is both. moie. simple and pio >- 
ably more realistic than the commonly used isolated sphere idealization. Since planar fronts 
present less evaporating surface area per unit mass of cool gas compared than do spherical 
clouds, the predicted densities of such intermediate temperature diagnostic ion species as U 
VI are below other models, and in better agreement with observations. Furthermore, the time 
dependence of these fronts is characterized by a rather lengthy period in which important ion 
column densities { e.g ., N V, 0 VI) evolve very slowly, although this depends somewhat on 
the magnetic field structure. These fronts are a significant source of UV emission in galactic 
haloes. 

Balbus and B. Dalton (graduate student on the grant) have developed an analytic model 
for the evaporation of a cold cloud in a hot gas in which the thermal conductivity vanes 
smoothly from a classical diffusive form to a nondissipalive, flux-limited form. Ibis model 
will be of considerable interest because the mathematical form of the flux-hunter used is 


4 



quite standard in numerical work, but to date no benchmark analytic problem has been 
available for code testing. 

Balbus and Ilawley have found that astrophysical accretion disks are generally unstable 
to axisymmctric disturbances in the presence of a weak fsublhermal) magnetic fields. 1 he 
identification of this instability may solve the long standing puzzle over the nature of the 
physical origin of the turbulent viscosity invoked as a transport mechanism in disks, as well as 
offer a remarkably simple way to efficiently amplify magnetic fields in rotationally supported 
systems. 

Borkowski, Balbus, and C. Fristrom (post-doc on the grant) studied the evolution of 
plane-parallel magnetized thermal conduction fronts in the interstellar medium (ISM). Sep- 
arating the coronal (~ 10 6 K) ISM phase and interstellar clouds, these fronts have been 
thought to be the site of the intermediate-temperature regions whose presence was inferred 
from 0 VI absorption line studies. The authors followed the front evolution numerically, 
starting from the initial discontinuous temperature distribution between the hot and cold 
medium, and ending in the final cooling stage of the hot medium. Initially, electron thermal 
conduction transports energy from the hot to the cold medium, evaporating and heating 
the cold gas to temperatures slightly over half of the initial hot medium temperature. The 
front broadens continuously until radiative losses balance conductive healing in the evapo- 
rated gas. Subsequently, the evaporated gas cools down and a condensing front propagates 
into the hot medium. For the typical ISM pressure of 4 x 10 3 Kcm -3 and the hot medium 
temperature of 10 fi K, the transition from evaporation to condensation in a nomnagnetized 
front occurs when the front thickness is 15 pc. This thickness is by a factor of 5 smaller than 
previously estimated. The 0 VI column densities in both evaporative and condensation 
stages agree with observations if the initial hot medium temperature 1\ exceeds 7.5 x 10* K. 
Condensing conduction fronts give better agreement with observed 0 VI line profiles be- 
cause of lower gas temperatures. Collisionally excited C and O lines in both evaporative and 
condensing fronts, as well as He II .‘301 A line in the evaporative front, could dominate the 
diffuse UV and EUV line emission. A fraction of the observed soft X-ray background may 
also be produced in a conduction interface. An ordered magnetic field reduces conductive 
energy transport, but the reduction is substantial only for large inclinations of magnetic field 
lines in the hot gas to the front normal. 

Borkowski, Sarazin, and Soker considered the interaction of a moving planetary nebula 
(PN) with the interstellar medium. The PN shell is compressed first in the direction of the 
stellar motion. This produces a di pole asymmetry in the surface brightness of the nebula, 
typically at a nebular density of ~ 40 cm -3 if the nebula is located in the galactic plane. In 
the later stages of the interaction, this part of the shell is significantly decelerated with respect 
to the central star, and the PN becomes strongly asymmetric in shape. T his distortion and 
the subsequent stripping of the nebular gas away from the central star typically occurs at a 
low nebular density of ~ G cm -3 . Borkowski, Sarazin, and Soker examined the morphology 
of PNs with central stars whose proper motions exceed 0.015 arcsec yr -1 , and found that 
many of the extended nebulae are interacting with the interstellar medium (ISM). Their 
sample doubles the number of known PNs interacting with the ISM. They also examined 
the morphology of nearby PNs, and found a number of strongl}' asymmet ric nebulae. I he 
expectation is that the central stars in these PNs should have large proper motions, and the 
authors predicted the directions from the observed asymmetries. They suggested the use of 
interacting PNs to determine the filling factor of the warm neutral and ionized components 
of the ISM. 

Brcgman and Houck completed their investigation of the behavior of low temperature 



galactic fountains. Observations show that neutral gas is present in the halo of the galaxy 
with a characteristic height of 1 kpc, Nj{ « 10 19 — 10 2 crn~ , and -100 kin/s < v < 100 km/s. 
To examine whether these observations can be reproduced by the galactic fountain model, 
analytic as well as one and two dimensional hydrodynamic calculations were performed. 1 he 
temperature and density at the base of the halo are the important free parameters. Fluid 
flow in the halo was calculated in a Cartesian frame with a constant gravitational force and 
optically thin radiative losses. When the cooling time t coo i is less than the sound crossing 
time /, in the halo gas, a transonic fountain occurs in which cool clouds are produced in the 
upward flow. If t, < t coo i < 10/,, a steady subsonic fountain occurs, also with a well-defined 
cloud formation layer. In these cases, the density at the base of the flow determines the 
velocity and mass flux, while the choice of temperatures determines the height of the cloud 
formation layer. For halo pressures appropriate for the solar circle, transonic flow occurs for 
T < 3 x 10 5 I< and n < 1.5 x 10 -3 cin -3 . When t coo i > 10/, (which occurs when T > 10 6 K 
and n < 5 x 10 -4 cm -3 ), convective or Rayleigh-Taylor instabilities disrupt the steady flow, 
broadening and complicating the region of cloud formation. The model that best reproduces 
the observations is near the transonic regime and has a temperature and density at the base 
of the fountain of T = 3 x lO 5 K and u = 1 X 10 -3 cm -3 . 1 he cooling time for the hot gas is 
t CO ol — 3-3 x 10 7 yr and the time for newly formed clouds to return to the disk is '1.7 x 10 7 yr, 
neglecting drag forces. If extended uniformly over the galaxy using a radius of 1U kpc, the 
mass flux rate onto one side of the disk is 0.2 A/ 0 /yr and the energy radiated is 2.9 x It) 39 
ergs/s, mostly in the extreme ultraviolet region. 

P. Plait (a graduate student on the grant) and Sokcr studied the density and mass profile 
of the planetary nebula NGC 0826. They discussed the possible formation of the dense shell 
at the edge of the nebula as a result of interaction with the ISM. The formation of the 
spherical outer halo together with an elliptical inner region can be explained by a low mass 
companion to the progenitor of NGC 0826 that was evaporated or collided with the central 
star. 

Sarazin and C. Grancy (a graduate student on the grant) have calculated the emission 
from gas at temperatures of 7’ ~ It)*’ K due to optical coronal lines that result from fine 
structure transitions within the ground states of highly ionized atoms. These lines include the 
[Fe X] A6374 line, the [Fe XIV) A5J03 line, and many other lines from stages of ionization 
of S, Ca, Fe, and Ni. The calculations are valid in the limit of very low density, and 
apply to the diffuse gas found in supernova remnants, the interstellar medium in galaxies, 
and intergalactic gas. The calculations are done for two physical situations. First, the 
optical coronal line emission is calculated for a gas in ionization equilibrium for temperatures 
between 10 5 and 10 7 K. The results were given as a function of the density and of t he pressure 
of the gas. Second, the emission of a gas cooling under its own radiation is calculated. I he 
emission is calculated for gas cooling at either constant pressure (isobarically) or constant 
density (isochorically). The effects of nonequilibrium ionization are included, but Grancy 
and Sarazin find these to be insignificant. They also find that the amount of optical coronal 
emission produced by a given amount of cooling gas is largely unaffected by changes in 
element abundances (providing the proportions of heavy elements remain the same). 

Soker studied the formation mechanisms for the two optical bright knots in planetary 
nebulae, which are called “ansae”. Soker used a two-dimensional hydrodynamic code to 
study the shaping of young asymmetrical planetary nebulae. He suggested that the shape 
of a young asymmetrical planetary nebulae may be influenced by a close binary-star system 
located at its center. This binary system is a relic of the common envelope phase, presumably 
through which the asymmetrical planetary nebula evolved, lie assumed that for a short 
period of time, shortly after the cease of the slow wind and long before the fast wind becomes 


0 



effective, the binary system ejects a small amount of mass, mainly in the equatorial plane. 
He found that at late times the high density region has a “horseshoe” shape, as viewed in 
the symmetry plane. There is an instability in the maximum density region. 

In order to provide a theoretical framework for observations of interacting PNs, Soker, 
Borkowski, and Sarazin calculated hydrodynamical models of PNs interacting with the ISM. 
In most cases, their results are in fair agreement with those obtained with the help of a tliin- 
shell approximation, and agree with observational data. 1 his simple picture breaks down for 
a high-velocity PNe in a low density environment, where the shocked ISM cools down slowly. 
In such a case, a Rayleigh-Taylor instability develops, leading to the shell fragmentation. 
Soker, Borkowski, and Sarazin discuss in this context a recently discoved PN in the globular 
cluster M 22, which is being stripped away from the cluster by the ram pressure of the ISM. 
Similar considerations are of importance in elliptical galaxies where PNs are stripped away 
from their central stars by the hot ISM, whose presence is inferred from X-ray observations. 

Soker and M. Livio (Teclmion, Haifa) used two dimensional hydrodynamics to investigate 
the problem of shaping of planetary nebulae by the interacting winds model. They obtained 
the shock structure and demonstrated that interacting winds are capable of producing mor- 
phologies similar to the ones observed in planetary nebulae, provided that a density contrast 
exists between the equatorial and polar directions. They discussed mechanisms that can 
produce the required density contrast and showed that binary central stars provide the most 
natural mechanism, especially via common envelope evolution. It is not yet entirely clear 
whether stellar rotation is sufficient to generate the required contrast in the cases of single 
central stars. They finally pointed out that interacting winds may play an important role in 
SN1987A. 


Hot Gas in Elliptical Galaxies 

Bregman and J. Boisseau (an undergraduate astrophysics major at Virginia) investigated 
a particle acceleration model with application to the A(JN environment. '1 hey considered 
a time dependent model in which plasma enters a shock region and undergoes first order 
Fermi acceleration. The time evolution of the particle distribution and the associated syn- 
chrotron emission has been calculated for an energy dependent diffusion coefficient. \Y hen 
the emitting region is not spatially resolved, the spectrum is a power-law' with an exponen- 
tial turnover in frequency. With increasing time, the power-law and the turnover extend to 
higher frequencies. To test the model, observations at three or more frequencies should be 
made of time-resolved flux increases from the optically thin nonthermal emission region in 
variable active galactic nuclei. Existing data that can be used to lest the model are scarce 
but offer qualitative support. 

Bregman and M. Roberts (NRAO) studied the III properties of cooling flow elliptical 
galaxies. The Arecibo 305 m radio telescope was used to make 21 cm measurements of NCIC 
4406, an elliptical galaxy containing 1U ,0 jUq of hot diffuse gas. 1 he III decreases by a factor 
of 3 one beam width from the center (3.5 aremin, or 20.4 kpc for d = 20 Mpc) and by an 
order of magnitude two beam widths from the center. Excess III emission is seen in the 
X-ray tail, but not in dust-rich regions. The III measurements do not support the suggestion 
that the X-ray tail is hot gas being stripped from the galaxy. The velocities of the peaks 
of the emission differ by less than 20 kni/s for observations separated by moie than 40 kpc, 
indicating that the III has little rotation. Infall is also inconsistent with the line pi Miles, and 
we argue that random motion of III clouds is the probable origin of the III profile shapes. 
The properties of the III gas agree with expectations of the cooling flow model in which (old 
gas is produced in situ as the hot gas radiative!}' cools. 


Sarazin and G. Ashe (an undergraduate astronomy major supported during a summer on 
the NASA grant) calculated steady-state cooling flow models for the hot gas in normal ellip- 
tical galaxies, including the loss of X-ray emitting gas through inhomogeneous cooling. Since 
this cool gas may eventually form into stars, they refer to this process as “star formation”. 
They find that the star-forming cooling flow models have smaller sonic radii, smaller inflow 
rates in their central regions, lower densities, and higher temperatures than homogeneous 
non-star-forming models. The star-forming models have higher velocities in the outer re- 
gions and lower velocities in the inner regions. Most of these differences are pronounced only 
in the central regions of the galaxy. The star-forming models also have lower X-ray lumi- 
nosities (by a factor of at most 2.6), and this brings their luminosities into better agreement 
with the average observed X-ray luminosities of elliptical galaxies. However, there is large 
dispersion in the observed X-ray luminosities of ellipticals, and this cannot be explained by 
variations in the efficiency of star formation. The star-forming models have X-ray surface 
brightness profiles which are much less centrally peaked than the non-star forming models, 
and models with very efficient star formation are in reasonable agreement with the observed 
surface brightness profiles of ellipticals. The integrated X-ray spectra of star -forming and 
non-star-forming models are rather similar. On the other hand, the X-ray spectra of the 
outer regions of the star-forming models have much stronger soft X-ray line emission, and 
this could be used to test between these two classes of models. 

Sarazin and ILay White (University of Alabama) have studied the wide dispersion in 
the X-ray luminosities of elliptical galaxies of a given optical luminosity. It is known that 
the X-ray luminosities of elliptical galaxies are strongly correlated with their optical lumi- 
nosity, but there is a considerable scatter in this relationship. They have searched for a 
correlation between the X-ray luminosity and many intrinsic properties of elliptical galaxies, 
and they have found no strong correlations. Instead, they find that the X-ray luminosities 
correlate strongly with environment. The less luminous X-ray ellipticals are found in dense 
environments; that is, for a given optical luminosity, the X-ray luminosity of an elliptical 
anticorrelates with the presence of many neighboring galaxies. Within the Virgo cluster, 
they also find that the X-ray luminosity increase with distance from the cluster center. This 
suggests that the amount of hot gas in elliptical galaxies is controlled, at least in part, by 
the stripping of gas due to ram pressure ablation or collisional stripping. 


Ilot Gas and Galaxy Formation 

Chevalier and E. Lufkin (a graduate student on the grant) examined t he spherical collapse 
of a gas cloud initially in hydrostatic equilibrium using self-similar solutions. They assumed 
an initial power law density distribution and that some energy loss of cooling initiates the 
collapse. The continued collapse can occur on a hydrodynamic timescale (hydro-regulated) 
or be regulated by the cooling process (cooling-regulated). The hydrodynamic timescale is 
roughly equal to the free-fall time. Collapses under self-gravity are generally hydro-regulated, 
while collapses in an external gravitat ional field are generally cooling-regulated. The collapse 
of gas clouds to form galaxies could occur on a hydrodynamic timescale only if the gas 
dominated the gravitational field and the inflow formed a central dominant mass. 1 he 
likelihood of a collisionless component and the observation of extended mass distribut ions in 
galaxies implies that the collapse of a hot prologalactic cloud is cooling-regulated. Chevalier 
and Lufkin examined how our Galaxy might form from an initially smooth, hot cloud and 
showed that some nonlinear substructure must initially be present. 

S. Murray (a post-doc on the grant) and D. Lin (U.C. Santa Cruz) have examined the 
diffusion of metals within protoglobular clusters The observed homogeneity in [Fe/Ill among 
the stars within a typical globular cluster implies that the cluster stars were formed in well 
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mixed protoglobular cluster clouds over a timescale short compared with that required for 
gas to recycle and self-enrich. Previous work has indicated that the protoclustcrs can form 
from primordial gas on a timescale comparable to the crossing time of the clusteis in the 
galaxy and undergo star formation through fragmentation on a timescale comparable to both 
the internal dynamical timescale of the clusters and the stellar evolution timescale of upper 
main sequence stars. In their current this work, they examine various transpoit processes 
which may have aHected the chemical homogeneity in the protocluster clouds. 1 hey find that 
the characteristic diffusion timescale associated with collisions between grains and gas atoms 
is considerably longer than that on which star formation is expected to occur. Collisions 
between large grains and gas atoms lead to mass segregation and metallicity gradients on 
a timescale comparable to the crossing time of the clusters in the Galaxy. One possible 
mechanism for inducing and maintaining chemical homogeneity is turbulent diffusion in 
the clouds. The mixing timescale required in this case is comparable to several internal 
dynamical timescales, longer than the evolutionary timescale of the most massive stars, ami 
shorter than the Galactic orbital timescale of the clouds. Thus, metals in presently observed 
stars probably did not originate from upper main sequence stars of a coeval generation. 


Hot Gas in Clusters of Galaxies 

Balbus studied the effects of a magnetic Held on the development of linear thermal 
instability in a hot plasma. He found that local condensational growth is possible for a 
class of short wavelength modes due to magnetic confinement of higher pressure regions. 
In regions of ordered magnetic Helds within clusters, this effect may be of (destabilizing) 
importance if thermal conduction can be suppressed. 

Balbus, Hawley and Eric Lufkin (graduate student) continued their work on the forma- 
tion and evolution of finite-amplitude and finite-wavelength disturbances in cluster cooling 
flows perturbed by a single galaxy orbiting along the symmetry axis. They have performed 
two-dimensional time-dependent axisymmelric numerical simulations and found t hat for den- 
sity and temperature profiles typical of a large cluster cooling flow, the perturbing galaxy 
excites nonlinear gravity modes which result in local density enhancements of order unity 
and fluid velocities iu excess of 500 km s _1 . 

Balbus and Murray (postdoc on the grant., now at Cambridge) completed a study of 
time dependent cluster cooling flows using a characteristic based scheme. I hey found that 
evolutionary efTccts may be important, because the effects lend to reproduce a linear mass 
accretion profile unrelated to any in situ cooling or deposition process. 

Balbus and Soker have concluded a study of linear thermal instability in optically thin 
dynamical and static plasmas. Several new and fundamental results have been found. For 
example, they have shown that systems in both hydrostatic and thermal equilibi ium aie 
in general thermally stable if convcctively stable, and vice-versa. Ibis basic ami important 
result went undiscovered for some 23 years since the original paper on thermal instability. 
Balbus and Soker have also examined dynamical flows in some detail. This work has been 
applied to the how of hot x-ray emitting gas observed in the volume between the galaxies in 
richly populated clusters. It has been shown that these flows are considerably more stable 
than commonly thought. This has caused investigators to reconsider current theories of 
evolution of these hot plasmas. 

Balbus and Soker concluded a study of the resonant interaction between the cubital 
motion of a galaxy in a cluster of galaxies and the buoyant oscillations of the hot X-ray 
emitting gas bound to the cluster. Such interactions excite internal gravity waves iu the 
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gas, (Internal gravity waves are well-known in oceanographic and atmospheric studies.) 
These waves represent the fundamental dynamical response of the gas to external forces 
with characteristic time scales of order galactic orbital periods. The presence of these waves 
has many consequences. The energy carried by these waves can be a significant source of 
drag for an orbiting galaxy. Angular momentum transport by the waves couples galactic 
orbital motion and the dynamics of the hot gas. Gravity waves propagate toward the central 
regions of the gas, where they are likely to grow to nonlinear amplitudes, and are perhaps 
associated with optical filamentary structure often observed near the center of clusters. The 
group velocity associated with internal gravity wave propagation is an important dispersive 
mechanism that limits the ability of thermal instabilities to grow in a hot, cooling plasma. 

P. Becker (a post-doc on the grant) and Balbus have investigated relativistic particle 
transport in spherical background flows with momentum and spatially dependent diffusion 
coefficients. They have obtained results that describe the diffusion and first-order Fermi 
acceleration of relativistic particles in a gas of scatterers experiencing either inflow or outflow. 
Four qualitatively different regimes of behavior have been identified, depending on the sign of 
the velocity of the background plasma and the competition between advection and diffusion 
at large and small radii. Solutions have been obtained for the Green’s function corresponding 
to monoenergetic particle injection at a thin shell. Particles propagating away from this 
shell always gain energy if the scattering centers are inflowing; otherwise they lose energy. 
Results have been used to calculate the synchrotron spect ra produced by relativistic electrons 
accelerated in cooling flows in clusters of galaxies. Particle acceleration in converging cooling 
flows may provide a natural explanation for the extended (j> 100 kpc) radio halos observed 
around the central galaxy in a number of cooling-flow clusters, most notably Perseus. The 
results of the transport theory can also be applied to the broadening and redshifting of line 
profiles in supernova winds and to the production of asymmetric line profiles in quasars. 

Bregman and L. David (former NASA theory post-doc now at Ilarvard/CFA) applied 
numerical hydrodynamics to investigate the effects of heat conduction and magnetic fields 
on the growth of thermal instabilities in cooling flows. For magnetic field strengths of about 
1 fxG (similar to that in inlraclustcr gas), magnetic pressure forces can completely suppress 
shocks from forming in thermally unstable entropy perturbations with initial length scales 
as large as 20 kpc, even for initial amplitudes as great as 00%. The suppression of shock 
formation in cooling condensations significantly reduces the predicted luminosity of optical 
and ultraviolet emission lines produced as thermally unstable gas cools and decouples from 
a cooling flow. They find that perturbations with initial amplitudes of 50% and initial 
magnetic field strengths of 1 /*G cool to 10' 1 K on a time scale which is only 22% of the 
initial instantaneous cooling time. Nonlinear perturbations can condense out of cooling 
flows on a time scale substantially less than the time required for linear perturbations and 
produce significant mass deposition of cold gas while the accreting inlraclustcr gas is still at 
large radii. Our models with heat conduction indicate that in order for entropy perturbations 
to condense out of cooling flows, the magnetic fields must be sufficiently tangled to reduce 
the Spitzer conductivity by a factor of approximately 100. Observational consequences have 
been considered of condensations in cooling flows. It is found that Faraday rotation can 
be generated by cold, dense, high magnetic pressure clouds lying along the line of sight to 
cluster radio sources. 

Bregman and David also examined whether viable models could be found for duster 
cooling flows that avoided the prediction of the standard model that gas was cooling at 
enormous rates. The two most important observations that a model must reproduce is the 
observed temperature gradient and the lack of young stars that would be produced if cooled 
gas normally formed stars. In the first model considered, drag heating by galaxies reheats the 
cluster gas and also suppresses thermal instabilities in certain temperature regimes. Numer- 
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ical hydrodynamic simulations show that it is possible to produce the observed tempera . ure 
gradient while significantly reducing the cooling rate for the gas, but only lor a narrow range 
of the galaxy heating parameter. In another model, the gas cools and begins to Form stars, 
which lead to supernovae that reheat the cooled gas to the ambient cluster temperature. 1 Ins 
model is viable if only high-mass stars are formed and if the supernovae are subluimnous. 


Bregman, B. McNamara (a graduate student), and CTConncIl investigated the in hared 
properties of cooling flow systems and the behavior of dust in a hot environment. Ill AS strip- 
scans were used to study the infrared emission properties of 27 dominant galaxies that lie in 
the center of clusters rich in hot gas. In this sample, 30% of the objects were detected at 12/2 j 
microns while 46% were detected at 60/100 microns. The ratio of the infrared to optical flux 
densities is an order of magnitude greater than for elliptical and SO galaxies. I he presence 
of infrared emission is not obviously correlated with the redshift z, the optical magnitude 
my, the absolute magnitude My, the presence of young stellar populations, emission line 
gas, the X-ray luminosity, or the cooling rate of hot gas A/. The infrared luminosities, which 
are comparable to the X-ray and blue optical luminosities (1U 14 — 10 ergss ), imply dust 
masses of (0.3 - 3) x 10 7 (7’ rf „,,/30 K)~ 5 ' 9 Af 0 , and gas to dust ratios < m0{T du$t /W K) ' • 
The dust grains are probably heated by collisions with electrons ill the X-ray emitting gas, 
a process that also causes grain destruction. A possible supply of dust is stellar mass loss, 
although this may not be entirely adequate to balance the inferred destruction rale. 


McNamara and O’Connell completed their study of the spectra of cD galaxies in cooling 
flows. They present AA3400-5I00 spectrophotometry of the nuclei of 13 galaxies in cluster 
cooling flows with estimated cooling rates of M ~ 5-350 A/ 0 /yr. '1 hey find spectral anoma- 
lies, with respect to a control sample, in 8 of the 13 objects in the form of strong [O II) 
emission and/or excess flux in the 3500-4000 A continuum. The anomalies lend to be as- 
sociated with larger predicted cooling flow effects, though the relationship is not very tight. 
The continuum anomalies are apparently due to the presence of massive OB stars that have 
formed from the cooling flows, rather than nonlhermal radiation or inelalhcity ellccts. In 
no case, however, is the excess consistent with a 100% efficient conversion to stars with a 
normal IMF. 


McNamara, Bregman, and O’Connell searched for II I in 11 cluster cooling flows with 
the Arecibo and the Green Bank 300 ft. radio telescopes. They detected II I in absorption 
against the radio continuum sources near the centrally-dominant galaxies in the clusters 
2A0335+096 and MKW3s. The II I absorption features are redshiflcd with respect to l he 
central galaxies, indicating that the II I is probably falling toward the central galaxies. II 
was not detected in emission. Upper limits to the II I-lo-dust ratio in four objects are found 
to be consistent with the Galactic mean gas-to-dust ratio. These II 1 detections, and one lor 
NGC 1275 reported by other groups, are the most direct evidence to date that the gas in 
cooling flows reaches the central galaxy. 

McNamara and O’Connell completed two analyses of the properties of galaxies in the 
center of massive cooling flows in clusters of galaxies. X-ray observations imply that mass 
drop-out from the flows occurs at the rate of several hundred solar masses per year m 
these galaxies. The only plausible sink for this inflow is star formation. I here is strong 
corroborative optical evidence for the existence of the cooling flows, but in no case does the 
amount of observed star formation correspond to more than about 30 /o of the predicted 
rate if the initial mass function is normal. The only viable repository for the inflowing mass 
therefore appears to be in the form of very low mass stars or smaller objects— t.e. brown 
dwarfs” or “jupiters”. O’CounclI and McNamara examined in some detail the prospects lor 
direct detection of such objects and concluded that they are not good with current optical 


or infrared instrumentation. 


Sarazin has shown that X-ray absorption lines may be detectable from clusters of galaxies, 
if there is a background quasar projected behind the cluster. The detection of such a resonant 
absorption line and of the emission from the cluster in the same line would allow the distance 
to the cluster to be determined, independent of the standard extragalactic distance scale. 
The basics idea is simple; the absorption from the gas is independent of distance as long 
as the cluster is in front of the quasar. The flux of X-ray line emission depends inversely 
on the square of the distance, so a comparison of the absorption and emission gives an 
estimate of the distance. This technique is very similar to the distance detei ruination based 
on the Zeldovich-Sunyaev microwave diminution. Sarazin found that the required X-ray 
observations should be possible with AXAF, but that it will be difficult to attain the level 
of accuracy required to provide a definitive determination of the Hubble constant. 


Sarazin and Graney have calculated the predicted luminosity and the surface brightness 
profiles of the optical coronal emission lines for models of cooling flows in galaxy clusters and 
in individual elliptical galaxies. The surface brightness profiles of the lines directly indicate 
where gas is cooling below X-ray emitting temperatures. In homogeneous cooling (low models 
in which all of the gas flows into the central regions, the coronal lines are very centrally 
concentrated. For inhomogeneous models in which gas cools out of the flow at large radii, 
the surface brightness profiles of the coronal lines are more extended, and the central s in face 
brightness is much lower. Existing observations apparently rule out homogeneous cooling 
flow models for several clusters if the cooling rates inferred from X-ray observations aie 
correct. If gas is indeed cooling at these rates, it should be possible to detect optical coional 
emission lines in many cooling flow clusters. If the lines are not detected, it would establish 
unambiguously that gas is not cooling at the rates inferred from the X-ray obset vations. If 
the lines are detected, their spatial distribution will show directly where the cooling gas is 
being deposited. 


Sarazin and R. White (Univ. of Alabama) have compared the X-ray images of M87- 
Virgo and NGC1275/Perseus clusters to models for cluster cooling flows which include star- 
formation. X-ray observations indicate that M87/Virgo and NGC1275/Pcrseus have cooling 
flows which arc associated with accretion rates of ~ 20 — 30 and ~ 300 — 500 A /0 yr , 
respectively. By deconvolving X-ray surface brightness profiles of these cooling flow dusters, 
a number of authors have found the accretion rates to decrease inwards, suggesting that 
the bulk of the mass flux is dropping out (presumably to form stars). However, cooling 
flow models which assume the mass (lux is constant with radius (so there is no star for- 
mation) have never been calculated self-consistently for these objects. Sarazin and W hite 
determined whether star formation is necessarily occurring in these cooling flows by calcu- 
lating constant-inass- flux models for all reasonable parameter space. They assumed that 
the inflow was steady-state, and that collisional heating and thermal conduction were in- 
significant. They found no constant-mass-flux models which were consistent with all of the 
relevant observations, so concluded that mass is indeed dropping out of these cooling (lows. 
They then considered a variety of star-forming models in which the X-ray emission due to 
cooling condensates is taken into account. For M87/Virgo they could find models in which 
the mass flux vanishes by the center. In NGC1275/Pcrseus, they find that the mass flux 
diminishes by £ 40% by the time the flow reaches - 20 kpc (the effective inner limit of the 

X-ray surface brightness data). The remaining mass flux within 20 kpc was - 300 A / 0 yr" , 
which is unacceptably high for the flow temperatures we derive. 


Sarazin and Wise have calculated the effects of electron scattering on cluster cooling 
flows. Cluster cooling flows have electron scattering optical depths which are estimated to 
be about 10 -2 . Many of the cluster dominant galaxies centered in cluster cooling flows have 
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active nuclei, and about 1% of the luminosity from the active nuclei should appear as diffuse, 
scattered radiation. Due to significant amounts of diffuse emission in other wavebands, this 
scattered radiation is most likely to be observable at radio wavelengths. Distinguishing 
between scattered radiation and other sources of diffuse emission poses a serious problem. 
However, the scattered emission should have the same spectral index as the central source, 
should be axisymmetric, and has a surface brightness distribution which can be predicted 
from the observed X-ray surface brightness. They show that detection of this scattered 
radiation can provide a measure of the distance to the cluster that is independent of the 
Hubble constant. Observations of the scattered surface brightness profile can also be used to 
test for dumpiness in the gas, for variations in the central radio source on time scales 10 
years, and for the degree of beaming of radiation from the central source. 

Saraziu and Wise also studied the relationship between 13 L Lacertae (BL Lac) objects 
and Fanaroff-Riley type I (FR I) radio galaxies. They showed that the strong beams of 
radiation present in BL Lac s anti FR Fs will be scattered and made visible by cooling (low 
gas. 

Soker, Bregman, and Sarazin have shown that the nonlinear perturbations which lead 
to thermal instabilities in cooling flows might start as blobs of interstellar gas which are 
stripped out of cluster galaxies. Assuming that most of the gas produced by stellar. mass 
loss in cluster galaxies is stripped from the galaxies, the total rate of such stripping is 
roughly A f[SM ™ 100A/Qyr which is similar to the rates of cooling in cluster cooling 
flows. It is possible that a substantial portion of the cooling gas originates as blobs of 
interstellar gas stripped from galaxies. The magnetic fields within and outside of the low 
entropy perturbations may help to maintain their identities, by. suppressing both thermal 
conduction and Kelvin— Helmholtz instabilities. Inhomogeneities in the iutracluster medium 
caused by the stripping of interstellar gas from galaxies can have a number of other effects 
on clusters. These density fluctuations may disrupt the propagation of radio jets through 
the intracluster gas, which may be one mechanism for producing Wide-Angle-Tail radio 
galaxies. 

Soker and Sarazin have studied the effects which iutracluster magnetic field can have on 
cooling flows in clusters. Even a very weak magnetic field, if it is tangled on small scales, 
can effectively inhibit transport processes in the gas, such as thermal conduction, viscosity, 
and the settling of heavy ions. Estimates and limits on the strength of the magnetic Helds in 
clusters suggest that the field is not dynamically important there. However, when gas cools 
and flows into the center of cooling flow clusters, the magnetic pressure. is greatly increased 
by the compression and shear in the flow. they show that the magnetic field may become 
very important in the inner regions of cluster cooling flows. Ihe magnetic pressure becomes 
comparable to the thermal pressure within a radius of typically 1 — 10 kpc. W ithin this 
region, reconnection of the magnetic field becomes efficient and magnetic energy is released 
at a rate of 10 41 - 10 43 erg s~\ Such models predict large values for the Faraday rotation 
of central radio sources in clusters. 
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